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1 Introduction 



In recent years, Generalized Parton Distributions have become one of the main topics of interest in 
hadronic physics (for recent reviews, see, e.g., [||, 0, 0, |, |^). GPDs are a natural bridge between 
exclusive processes, such as elastic scattering, described in terms of form factors, and inclusive ones, 
described in terms of structure functions. As it happens for the usual Parton Distributions (PDFs), 
the measurement of GPDs allows important tests of non-perturbative and perturbative aspects of the 
theory, QCD, and of phenomenological models of hadrons. Besides, GPDs provide us with a unique 
way to access several crucial features of the structure of the nucleon. In particular, as pointed out 
first by Ji @, , by measuring GPDs a test of the Angular Momentum Sum Rule of the proton 
could be achieved for the first time, determining the quark orbital angular momentum contribution 
to the proton spin. Therefore, relevant experimental efforts to measure GPDs, by means of exclusive 
electron Deep Inelastic Scattering (DIS) off the proton, are likely to take place in the next few years 
P, [TO]. Besides, GPDs measurements will be done soon, also through fitting to the available data 



of the HI and ZEUS Collaborations, with the HERMES results used as cross checks but excluded 
from the fit. In this scenario, it becomes urgent to produce theoretical predictions for the behavior 
of these quantities. Several calculations have been already performed by using different descriptions 
of hadron structure: bag models [0, |12|, soliton models P, [l^, light-front approaches and 



phenomenological estimates based on parametrizations of PDFs |15, 16]. Besides, an impressive 
effort has been devoted to study the perturbative QCD evolution ]|l^ |18| of GPDs, and the GPDs 
at twist three accuracy [0]. 

So far, to our knowledge, no calculations have been performed in a Constituent Quark Model 
(CQM), although a step towards this can be found in ]Q where the non relativistic hmit is 
shortly discussed. The CQM has a long story of successful predictions in low energy studies of the 
electromagnetic structure of the nucleon. In the high energy sector, in order to compare model 
predictions with data taken in DIS experiments, one has to evolve, according to perturbative QCD, 
the leading twist component of the physical structure functions obtained at the low momentum scale 
associated with the model. Such a procedure, already addressed in |22|, has proven successful 
in describing the gross features of standard PDFs by using different CQM (see, e.g., ]|23|). Similar 
expectations motivated the present study of GPDs. In this paper, a simple formalism is proposed to 
calculate the quark contribution to GPDs from any non relativistic or relativized model. By using 
such a procedure, the GPDs can be easily estimated, providing us with an important tool for the 
planning of future experiments. 

The paper is structured as follows. After the definition of the main quantities of interest, the 
proposed calculation scheme is introduced in the third section. Then, results obtained in a simple 
harmonic oscillator model and in the Isgur and Karl model are shown in the following section. 
NLO perturbative evolution from the scale of the model to the experimental one has been evaluated, 
and results are presented in the fifth section. Further applications of the procedure, using relativistic 
models and including antiquarks degrees of freedom, are in progress and will be presented elsewhere 
||25|| . Conclusions are drawn in the last section. 
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2 General formalism 



We adopt the formalism introduced by Ji, who called GPDs "Off-forward parton distributions" in 
1^. The connection of these quantities with the "Non diagonal" ones introduced in 0] is discussed 
in |l^ and can be easily obtained. 



We are interested in diffractive DIS processes. The absorption of a high-energy virtual photon by 
a quark in a hadron target is followed by the emission of a particle to be later detected; finally, the 
interacting quark is reabsorbed back into the recoiling hadron. If the emitted and detected particle 
is, for example, a real photon, the so called Deeply Virtual Compton Scattering process takes 
place. Let us think now to a nucleon target, with initial and final momenta P and P', respectively. 
The GPD H{x, ^, A^), the main quantity we deal with in the present paper, is introduced by defining 
the twist-two part of the light-cone correlation function 

dX An Xn 



H{x,^,A')U{P')YU{P) + E{x,^,A')U{P')^^UiP) + ... (1) 



where A = P'—P is the momentum transfer to the nucleon, ellipses denote higher-twist contributions, 
ip is a. quark field and M is the nucleon mass. In obtaining Eq. (1), a system of coordinates has 
been chosen where the photon 4-momentum, = (go, q), and P = {P + P'Y /2 are collinear along 
z. The ^ variable, the so called "skewedness" , parameterizing the asymmetry of the process, is 
defined by the relation ^ = —n ■ A, where n = (1, 0, 0, — 1)/(2A) and A depends on the reference 
frame, being, for example, A = M/2 in the nucleon rest frame. The ^ variable is bounded by and 
Jap - Ay A. Besides, one has t = A^ = Ag - A^. 



In the r.h.s of Eq. (1), the dependence of the light-cone correlation function on the GPDs 
H{x,^,A'^) and E{x,^, A"^) is explicitly shown. By replacing, in the above equation, j'^ with the 
proper Dirac operator, similar expressions can be derived for defining polarized or chiral odd GPDs 
0. In the following we will only discuss the unpolarized, chiral even, twist-two GPD A^). 
As explained in [0, 0, unlike the usual PDFs, which have the physical meaning of a momentum 
density in the Infinite Momentum Frame (IMF), GPDs have the meaning of a probability amplitude. 
They describe the amplitude for finding a quark with momentum fraction x + (in the IMF) 
in a nucleon with momentum (1 -|- ^/2)P and replacing it back into the nucleon with a momentum 
transfer A'^. Besides, when the quark longitudinal momentum fraction x of the average nucleon 
momentum P is less than — ^/2, GPDs describe antiquarks; when it is larger than ^/2, they describe 
quarks; when it is between — ^/2 and ^/2, they describe qq pairs. 

The region \x\ > ^/2 is often called the DGLAP region, since the evolution of GPDs is 
governed there by the DGLAP equations of perturbative QCD |2^; the region |x| < ^/2 is called 



the ERBL region, because there the evolution is ERBL-like |2^. One should keep in mind that, 
besides the variables x, ^ and A^ explicitely shown, GPDs depend, as the standard PDFs, on the 
momentum scale at which they are measured or calculated. For an easy presentation, such a 



'^Note that in going from Ref. to Ref. Ji has redefined by ^. 
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dependence will be omitted in the next two sections of the paper, while it will be properly discussed 
in the last one. 

There are two natural hmits for H{x,^,A'^) : i) when P' — P, i.e., = ^ = 0, the so called 
"forward" or "diagonal" hmit, one recovers the usual PDFs 

H{x,0,0) = q{x) ; (2) 

ii) the integration over x is independent on ^ and yields the Dirac Form Factor (FF) 

j dxH{x,^,A')^F,{A') . (3) 

Any model estimate of the GPDs has to respect these two crucial constraints. 

3 A non relativistic scheme 

Our aim now is to evaluate the Impulse Approximation (lA) expression for ^, A^), suitable to 
perform CQM calculations. 

Let us start from Eq. (1). Substituting the quark fields in the left-hand-side one has 

J ATT 

r,r',k 

Y[br{k)ur{k)e-'^-^ + d^{k)vr{k)e''-'r]\p) 
and, taking into account the quarks degrees of freedom only, it becomes 

C = V / ^e^^"(P>,.(fc + A)6+(^ + A)e-^('=+^>^7^6,(^)M,(fc)e-*'=-^|P) . 

-J ZTT 

r,r',k 

Using lA and integrating over A yields 



C - EE /^e^'^^-t-(2'=+^))(P>,,(^+A)6+,(^+A)7''fe,,,(^)ii,(^)|P) = 

* r,r',k 

= E E ^ ('^ - ^ ■ (2fc + A)) (P>,,(^+ A)6+ ,(^+ A)7'^6,,,(^)xi,(^)|P) 



* r,r',k 

Let us introduce ^ = —n • A, so that Eq. (1) reads now: 
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E E <^ + f - ^ ■ ^) {P'\Urik + A)b+Ak + ^)Yb,^r{k)Ur{k)\P) = 
* r,r',k ^ ' 

= H{x, A')U{P')YU{P) + E{x, A')U{P')^^U{P) , (4) 

which holds exactly if the antiquark degrees of freedom are not considered. In fact, the l.h.s. is 
evaluated in lA and the r.h.s. is the leading-twist part of the light-cone correlation function, so that 
they have the same physical content. 

By taking the zero-components in the left and right hand sides of Eq. @), considering a process 
with < M^, one immediately sees that the contribution of the term proportional to ^, A^), 
in the right hand side of Eq. (^, becomes negligibly small, so that H{x, ^, A^) is given by: 

Hix,^,A') = ^ ^ 5 + I - (P>+(^ + A)6+ + A)6,,(^)«,(fc) |P) , (5) 



r'r,k 



where = ko + k^ has been introduced. In order to evaluate this expression by means of a CQM, 
one has to relate it to nucleon wave functions. In a non relativistic framework, if the normalization 
of the nucleon states is chosen to be 

(P'|P) = (27r)=^<5(P' - P) , 

for a symmetric wave function (as is the case in a quark model once color has been taken into 
account), one has (see, e.g. , p8| ) 



EE(^>^(^ + ^)^tr'(^+^)kr(^)Mr(^)|P) = 3 J i^* {h, h, k + A)ij{ki, h, k)dkldk2 = 

= h{k,k + A) (6) 

where the one-body non diagonal charge density 

p{f,r') = ip*{ri,r2,r')i^iri,r2,r) dri dr2 (7) 



and the one-body non-diagonal momentum distribution n{k, k + A) have been introduced. In terms 
of the latter quantity, Eq. (|^) can be written 

H{x, ^,A') = Jdk 6(^x+^-- ^ j n{k, k + A). (8) 

The above equation, which is our basic result, allows the calculation of H{x,$,,A'^) in any CQM, 
and it naturally verifies the two crucial constraints, Eqs. (2) and (3). In fact, the unpolarized quark 
density, q{x), in the lA is recovered in the forward limit when A^ = ,^ = 0: 



q{x) = H{x, 0, 0) = J dknik) ^ ' J^] (9) 



so that the constraint Eq. (2) is fulfilled. In the above equation, n{k) is the momentum distribution 
of the quarks in the nucleon: 

n{k) = I d^^''-'"^ p{r,r')drdf^ . (10) 



As is well known, the relation between the quark momentum distribution and the quark unpolar- 
ized density, Eq. (P), can be found by analyzing, in lA, the handbag diagram, i.e., the leading twist 
part of the full DIS process (see, e.g., ||2^, |29|), assuming that the interacting quark is on-shell. So, 
from Eq. (|^), derived as the non relativistic reduction of the light-cone correlation function in the 
lA, the quark density appears as the forward limit. Besides, integrating Eq. (^ over x, one trivially 
obtains 



where p{r) = lim^^^p(?^, r) is the quark charge density. The r.h.s. of the above equation gives the 
lA definition of the charge FF 



dfe^^Xf) = F{A') , (11) 

so that, recalling that F(A^) coincides with the non relativistic limit of the Dirac FF Fi(A^), the 
constraint Eq. (3) is immediately fulfilled. 

Let us introduce now the following sets of conjugated intrinsic coordinates 



^= -^(n + ri + ri) ^ K = ^{ki + k2 + k3) 
I ,^ ^ - 1 



P=—{rl-r2) ^ kp = —{ki-k2) 



A = -^(rl + r^2 - 2rl) ^ kx = ^{k^ + k2 - 2k^) 



in terms of which, in a coordinate system where = 0, the FF Eq. ([TT]) can be written [^, Q 

F(A2) = j dkpdkxr (kp,kx - ^/273A) ^ (^p, ^1) (12) 
and iJ(x, ^, A^), Eq. (||), by substituting Eq. (|^ into Eq. (||), reads 

i/(x,e,A2) = j dkpdkxr{kp,kx-^?>A^ij{kp,kx)5l^ + ^-^'^ . (13) 

One immediately realizes that Eq. (0) is obtained from Eq. (|T^) by performing the x integration. 

With respect to Eq. (||), a few caveats are necessary. First of all, due to the use of CQM 
wave functions, which contain only constituent quarks (and also antiquarks in the case of mesons), 
only the quark (and antiquark) contribution to the GPDs can be evaluated, i.e., only the region 
X > E,/2 (and also x < — ^/2 for mesons) can be explored. In order to introduce the study of the 
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sea region (— < x < ^/2), the model has to be enriched. To this respect, calculations including 



a substructure of the constituent quark, as proposed by several authors ||31], |32|, |3^, are in progress 
and will be presented elsewhere p5|. 



Secondly, we remind that Eq. holds under the condition ^ M^. If one wants to treat 
more general processes, such a condition can be easily relaxed by keeping the terms of 0(A^/M^) 
in going from Eq. @) to Eq. (||). At the same time, an expression to evaluate ^, A^) could be 
readily obtained. 

Finally, in the argument of the 5 function in Eq. (10), due to the approximations used, the x 
variable is not defined in its natural support, i.e. it can be larger than 1 and smaller than ^/2. Several 
prescriptions have been proposed in the past to overcome such a difficulty in the standard PDFs case 
[p2| , |23[] . The support violation is small for the calculations that will be shown here. However one 
has to be cautious in interpreting the results after pQCD evolution is perfomed. A deeper discussion 
of this issue is beyond the scope of the present work P^ . 

We stress that our definition of if(a;, A^) in terms of CQM momentum space wave functions 
can be easily generalized to other GPDs, and the relation of the latter quantities with other FFs (for 
example the magnetic one) and other PDFs (for example the polarized quark density) can be 
recovered. Therefore the proposed scheme allows one to calculate GPDs by using any non relativistic 
or relativized |Q CQM, and it is also suitable to be implemented by corrections due to a possible 
finite size and complex structure of the constituent quarks, as proposed by several authors ^ ^]. 



4 Results in non relativistic quark models 

As an illustration, in this section we present the results of our approach in the CQM of Isgur and 
Karl (IK) In this model the proton wave function is obtained in a one gluon exchange potential 
added to a confining harmonic oscillator (h.o.) one; including contributions up to the 2huj shell, the 
proton state is given by the following admixture of states 

|A^) = as\^Sii2)s + 0'S'?S[i2)s + a^p5'i/2)M + aB|^£'i/2)M , (14) 

where the spectroscopic notation with t = A,M,S being the symmetry type, has been 

used. The coefficients were determined by spectroscopic properties to be |3^: as = 0.931, ag' = 
-0.274, ttM = -0.233, ax, = -0.067. 

If = 1 and ag/ = = a© = 0, the simple h.o. model is recovered. Let us now calculate the 
GPD H in the IK model by using Eq. (0). The different components appearing in the momentum 
space wave functions, obtained from Eq. (|1^ in the IK model, can be found in for the h.o. 



model, the corresponding wave function in momentum space reduces to p9| , ^ | 



ilj{kp,kx) = o/^ o , (15) 



„2 _ 1 QC -f-2 ; 



where the h.o. parameter can be fixed to a = 1.35/^ in order to reproduce the low t behavior of 
the charge FF, i.e., the r.m.s. value of the proton radius. 

The results in the IK model for the GPD H{x, ^, A^), for the fiavours u and d, respectively, neglecting 
in (inD the small D-wave contribution, are found to be: 
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M ( 3 



27ry 



e 3c«2 



cZ/Ca; / dky fo^kx, ky, X, ^, A. ) 



X 



(^fs{kx, ky, x,^, A ) + /(/cx) ^2/) X, ^, A ) j 



(16) 




fo{kx, ky, X, ^, A^) 



kp 
ko + k; 



'fa{Ax,kx)fa{Ay,ky)fa{Az,kz) , 



kz — 



M\x + ^/2f-{m^ + kl + ey) 
2M(x + ^/2) 



(18) 
(19) 
(20) 



fsikxi ky, X, ^, A 



' 2 2 

-a, + a„, 



6 q;2 + 2 q;4 + 3^2 1 3 +^-^'^) I ^2 



- 1 



+ a' 



M 



5 1A;2 lA;^ 2 k 9^^ - 



3q;2 \^ 3 
2 



+ asas' 



1 - 



k' 



2 /A^ 



3q;2 I 3 



A- A; 



(21) 



f{kx,ky,x,^,A'^) = -0505/-^ 



3a2 I 3 



— aMdS' 



+ 



V2 
2^ /A2 



6 a2 + 2 ^4 3a2 1 3 



3^4 V 3 



— + A-k 



(22) 



and kp — yjrn^ + + + the interacting quark has been assumed to be on shell. 

A few comments are in order: 
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The x-integration of Eqs. (|T^) and ([T7|) give the u and d flavor contribution to the proton 
charge FF in the lA, respectively, as given in 



In the forward limit, ^ = 0, = 0, Eqs (0) and ([T7|) give the distributions Hu{x, 0, 0) = u{x) 
and Hfi{x, 0, 0) = d{x) in the Isgur and Karl model, according to the findings of Ref. 



The results in the simple h.o. model can be immediately found from the ones presented above, 
just putting as = 1, as' = a_M = in Eqs. (^) and (|2^). In particular, by using the wave function 
Eq. ( |l5l) in Eq. (|l2l), one gets trivially 



(23) 



Besides, taking the "forward" limit, = ,^ = 0, of Eq. ([T3|) and substituting in the w.f. Eq. ([151), 
or, which is the same, taking the forward limit of Eqs. (l2l|) and ( p2t ) with as = 1, as' = o-m = 0, 
and performing analitically the integrations in Eqs. (16) and (0) one easily obtains: 



, 271 M / 3 \ 

Hix,o,o) = —r{-J 



dkke 2q 



k-(x) 



(24) 



where the integration limit k-{x) is 



k_(x) 



M 

T 



X 



ml 1 



M2 



X 



(25) 



and rriq is the quark mass. This is the same expression obtained in ||2^ starting from Eq. (9), using 



n{k) corresponding to the present model (called "model 1" in |^ 



n{k) 



_ 3k-' 
2a^ 



a- 



(26) 



Results are shown in Figs. 1 to 5. 

The behavior of the proton charge FF in the IK model is shown in Fig. 1. It is known |3^ that 
such a FF underestimates the data for large values of A^. 

Results for the GPD H at the low momentum scale fi^ associated with the model, are shown in 
Figs. 2 to 5. A value of rriq ~ M/3 has been used for all the estimates. Note that all the results 
are shown in the x > C,/2 region. We reiterate in fact that what has been calculated so far is the 
valence quark contribution to the full GPD H, so that we can provide estimates only in the positive 
DGLAP region. 

In Fig. 2, the A^ dependence of our results for the u and d flavors is shown. The forward, 
A^ = ^ = limit corresponds to the full line. One immediately realizes that a strong A^ dependence 
is found, in comparison with other estimates, for example with the one predicted in |jll[. This has 
to do with the already discussed strong t dependence of the FF in the IK model. 

In Figs. 3 and 4 we have the full A^ and ^ dependences. These findings, particularly clear 



from the three-dimensional plots in Fig. 4, 
dependence is a little stronger. 



are similar to the ones obtained in |]TT|, although the ^ 
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The scale of the left and right panels in Figs. 2 and 3 has been chosen in such a way that they 
would look exactly the same if the model were SU(6) symmetric (in fact, in that case, one would have 
Hu = 2ifrf). The observed difference clearly shows to what extent the 5*^7(6) symmetry is broken in 
the IK model. 

In Fig. 5, the comparison between the predictions of the IK and of the simple h.o. models is 
shown. As an example, results are presented for the d distribution at ^ = 0. 



5 QCD evolution of the model calculations 

According to a well estblished |22[, widely used scheme (see, for example, [Q), the results shown 
so far for H{x,^,A^) correspond to the low momentum scale /i^ associated with the model, and in 
order to compare them with the data of future experiments, one has to evolve them to experimental, 
high-momentum scales. We next proceed to do so. 

As already mentioned in the Introduction, an impressive effort has been devoted to studies of 
QCD evolution properties of GPDs, an essential feature to understand their physical content and to 
obtain the correct information from experiments. The QCD evolution of GPDs is presently known 
in both the DGLAP and ERBL regions, up to NLO accuracy [1T7|, 



The evolution of the results presented in the previous section has been carried on by using an 
evolution code kindly provided by Freund and Mc Dermott, adapted by us to our specific case. The 
features and performances of such a code are described in |]l8l, and an interface to it is available at 



the web site In order to be used as input in the evolution code, our GPDs have been translated 



into the Golec-Biernat, Martin off-diagonal parton distributions jF(xi,C) El- The new distribution 
and the variables Xi and ( are obtained from our quantities H, x and ^ according to definitions given 
in . Once the evolution has been performed, results have been translated back to our notation to 
allow a consistent presentation. 

The scale to be associated with the model is a low one, and the choice of its value is part of the 
model assumptions. We have chosen here /i^ = 0.346*6^^^, corresponding to the initial scale of the 



so called valence scenario of ||38|, since our input at the scale of the model is given by the valence 
quark contribution only. A thorough discussion about the choice of the initial scale can be found in 
[p3| . Of course, being the starting scale so low, it becomes very important to perform the evolution 
as accurate as possible. From this point of view, the NLO level, used here, allows a safe result. 

The evolved GPDs at = 10 GeV^ are shown in Figs. 6 and 7, for two values of ^ and A^, for 
the Non Singlet (NS) part of the u and d flavor distributions, together with the input at the initial 
scale. One should notice again that only the positive DGLAP region is shown, for both the input 
and the evolved distributions. As happens in the conventional PD case, and is easily understood, 
evolution lowers the mean x of the distribution, i.e., the partons accumulate near x = ^/2. 

In Fig. 8 a three-dimensional plot of the evolution, for fixed = 0.1 and fixed = —0.2 
GeV^, is also shown. 
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6 Conclusions 



Generalized Parton Distributions (GPDs) are a useful tool to access several relevant features of 
the structure of the nucleon, such as the angular momentum sum rule. A systematic theoretical 
study of many aspects of these objects started few years ago and it is being carried on in these 
years. The future experimental effort to measure GPDs is also ambitious, and, in this respect, 
theoretical estimates will be necessary for the planning of future experiments. In the present paper 
we propose a general formalism to investigate GPDs by means of non relativistic or relativized 
Constituent Quark Models. Starting from the general field-theoretical definition of the related light- 
cone correlation function, by performing an Impulse Approximation analysis and the non relativistic 
limit, the unpolarized, leading-twist GPD ^,A^) is obtained in terms of the nucleon wave 
functions in momentum space. From its expression, the quark momentum density is recovered as the 
forward limit, and the charge Form Factor as its x integral. Results for the valence quark contribution 
to ^, A^), in a simple harmonic oscillator model, as well as in the Isgur and Karl constituent 
quark model, are shown to have the general behavior obtained in previous estimates. NLO evolution 
to high experimental scales of the low momentum results obtained in the model has been performed. 
A proper treatment of the ERBL region within a constituent picture is presently under investigation. 

The proposed approach can have many interesting developments, such as the calculation of other 
GPD functions and DVCS observables, the use of relativistic models and the addition of effects due 
to a possible finite size and complex structure of the constituent quarks, as proposed by several 
authors. 
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Captions 



Figure 1: The Proton Charge Form Factor, as given by Eq. ( p!2D in the IK model. 

Figure 2: The GPD Hq{x,^, A"^), Eqs. ([T6| ) and ([T7|) , calculated for ^ = and three values of 
A^: the full line corresponds to = GeV^, the dotted one to A^ = —0.2 GeV^ and the dashed 
one to A^ = —0.5 GeV^. Left(right) panel: the u (d) flavor distribution. Notice that, due to the 
chosen scale, the left and right panels should look the same if the model under scrutiny were SU(6) 
symmetric. 

Figure 3: The GPD H{x,^,t), Eq. (21), calculated for A^ = -0.2 GeV^ and three values of 
^: the full line corresponds to ,^ = 0, the dotted one to ^ = 0.1 and the dashed one to ^ = 0.2. Left 
(right) panel: the u (d) flavor distribution. Notice that only the x > ^/2 region is shown. Due to 
scale used, the left and right panels should look the same if the model under scrutiny were SU(6) 
symmetric. 

Figure 4: The GPD H{x,^,t), Eq. (21), calculated at fixed ^ = 0.1 and -A^ in the range — 
GeV^ - 0.5 GeV^ (left panel), and at fixed A^ = -0.3 GeV^ and ^ in the range - 0.3 (right panel). 



Figure 5: The GPD Hd{x,0, A'^), calculated within the IK (full) and the h.o. (dashed) models, 
for A^ = GeV^ (top), -0.2 GeV^, -0.5 GeV^ (bottom). 

Figure 6: The NS, q flavor GPD, H^^x, ^ = 0.1, A^ = -0.2 GeV^, = 10 GeV^) (full), evaluated 
evolving at NLO the initial IK model distribution H^^{x,^ = 0, A^ = -0.2 GeV^, = 0.34 GeV^) 
(dotted). Left (right) panel: the u (d) flavor distribution. Only the x > region is shown. Due 
to the chosen scale, the left and right panels should look the same if the model under scrutiny were 
SU(6) symmetric. 

Figure 7: As in Fig. 6, but at ^ = 0.2, A^ = -0.5 GeV^. 

Figure 8: The NS, u flavor GPD, H^^{x,i = 0.1, A^ = -0.2 GeV^Q2), with the range 

= 0.34 GeV2 - = 10 GeV^. 
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